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ABSTRACT 

Black hole mass determination in active galaxies is a key issue in understanding various lumi- 
nosity states. In the present paper we try to generalise the mass determination method based on 
the X-ray excess variance, successfully used for typical broad line Seyfert 1 galaxies (BLS1) 
to Narrow Line Seyfert 1 (NLS1) galaxies. NLS1 galaxies differ from BLS1 with respect to 
several properties. They are generally more variable in 2-10 keV energy band so the natural 
expectation is the need to use a different scaling coefficient between the mass and the variance 
in these two types of sources. However, we find that such a simple approach is not enough. Al- 
though for majority of the 21 NLS1 galaxies in our sample a single scaling coefficient (larger 
by a factor 20) provided us with a satisfactory method of mass determination, in a small sub- 
set of NLS 1 galaxies this approach failed. Variability of those objects appeared to be at the 
intermediate level between NLS1 and BLS1 galaxies. These exceptional NLS1 galaxies have 
much harder soft X-ray spectra than majority of NLS1 galaxies. We thus postulate that the 
division of Seyfert 1 galaxies into BLS 1 and NLS 1 according to the widths of the H/3 line is 
less generic than according to the soft X-ray slope. 
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1 INTRODUCTION 

Active Galactic Nuclei (AGN) are strongly variable in X-ray band. 
In a fixed time-scale quasars containing very massive black holes 
vary less than much less massive Seyfert galaxies so the X-ray vari- 
ability was used to estimate the black hole mass (e.g. Hayashida et 
al. 1998; Lu & Yu 2001, Markowitz & Edelson 2004, Papadakis 
2004). However, if we aim to have a mass determination method 
which works relatively accurately we must take into account that 
the black hole mass is not the only parameter of the accreting sys- 
tem. The second, important parameter is the dimensionless accre- 
tion rate which may affect the variability. 

In the literature, the three possibilities were suggested for the 
exact scaling of the X-ray variability of an AGN with mass and the 
accretion rate. 

The first and the oldest one is the claim that the variability 
scales only with the X-ray luminosity (Barr & Mushotzky 1986, 
Lawrence & Papadakis 1993, Green, McHardy & Lehto 1993, Nan- 
dra et al. 1997, Leighly 1999a, Markowitz & Edelson 2004). It was 
recently argued by Liu & Zhang (2008) at the basis of a sample of 
14 objects with black hole mass obtained directly from reverbera- 
tion that the mass itself does not affect the variations. 

The second possibility is a continuous two-parameter relation 
between the X-ray variability and the black hole mass and the bolo- 
metric (not X-ray !) luminosity found by McHardy et al. (2006). 
About this idea see also papers of Markowitz et al. 2003, Bian & 
Zhao 2003, McHardy et al. 2005, O'Neill et al. 2005. 



The third possibility is that there is no continuous scaling with 
accretion rate but instead there are two separate scaling laws with 
mass: one for Narrow Line Seyfert 1 (NLS1) galaxies and another 
one for Broad Line Seyfert 1 (BLS1) galaxies (Hayashida et al. 
1998, Turner et al. 1999, Lu & Yu 2001, Uttley, McHardy & Pa- 
padakis 2002, Nikolajuk, Papadakis & Czerny 2004). This possi- 
bility is an interesting option since the two classes of objects differ 
considerably. 

The Narrow Line Seyfert 1 (NLS1) galaxies were introduced 
as a separate class by Osterbrock & Pogge (1985). These objects 
had all properties in the optical band of the Seyfert 1 galaxies but 
their line widths were surprisingly narrow in comparison to typ- 
ical Seyfert Is. The original division of Seyfert galaxies, accord- 
ing to the width, into narrow and broad line objects (Seyfert 1 
and Seyfert 2) seemed to be accompanied by the increase of the 
[OIITJA5007 to the H/3 flux ratios (Shuder & Osterbrock 1981). 
However, a class of objects with narrow lines but faint [O III] A5007 
was found and a new class had to be introduces. 

Till now the NLS1 as a class are defined by the requirement 
that the [OIII]A5007 to the H/3 flux ratios are smaller than 3 and 
Full Width Half Maximum (FWHM) of H/3 line in NLS Is was 
less than 2000 kms -1 . Further observations revealed other pecu- 
liar properties. Blue wings of some emission lines were stronger 
compared to that in Seyfert Is. Also the Fell lines are more asym- 
metric and much stronger than in Seyfert Is (Stephens 1989; Boro- 
son & Green 1992; Boiler et al. 1996; Wang et al. 1996; Leighly 
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1999b; Veron-Cetty & Veron 2000; Grupe et al. 2004). The Bald- 
win effect, which is common in BLSls, was seen in those objects 
as well (e.g. Warner et al. 2004; Shields 2007). However, Leighly 
& Moore (2004) observed the UV lines and found that the Equiva- 
lent Widths of emission lines in NLSls were systematically offset 
to lower values at a given continuum luminosity in comparison to 
BLS1 galaxies. 

Further differences came from X-ray observations. NLSls 
very frequently exhibit rapid and/or high-amplitude X-ray vari- 
ability. They systematically show larger amplitude variations than 
BLS1 on time-scales of minutes, hours and days (Hayashida et 
al. 1998; Fiore et al. 1998; Leighly 1999a; Turner et al. 1999; 
Markowitz & Edelson 2004; Uttley & McHardy 2005). Giant- 
amplitude X-ray variability has been observed in several NLSls 
(e.g. in NGC 4051 - Uttley et al. 1999; in PHL 1092 - Brandt 
et al. 1999). The X-ray spectra of NLSls are generally steeper 
than spectra of BLSls (Brandt, Mathur & Elvis 1997). The aver- 
age photon index, r2-iok c v in BLS1 in 2-10 keV band is ~ 1.7 
(Nandra & Pounds 1994; Reynolds 1997), whereas NLS1 galaxies 
have r 2 -iokov ^ 2.2 (Leighly 1999a). Similar effect is in 17-100 
keV band. The photons index in BLS1 galaxies is ri7_iookcV = 
2.02 ± 0.16 and in NLS1 V 17-100keV — 2.6 ± 0.3 (Malizia et 
al. 2008). The effect is much stronger in the soft X-ray band, 0.1- 
2.4 keV. The photon index, ro.1_2.4kcV, of BLSls varies around 
value 2.1-2.3 (Boiler et al. 1996; Grupe et al. 1998; Pfefferkom et 
al. 2001) but in NLSls it is typically of order of 3.0-3.3, and oc- 
casionally even larger (Boiler et al. 1996; Lawrence et al. 1997; 
Veron-Cetty et al. 2001; Williams et al. 2004; Bian 2005). More 
precisely, NLSls may show both steep as well as flat X-ray spec- 
tra, while BLSls always have flat spectra (Grupe et al. 1999). Gen- 
erally, there is an overall anticorrelation between the soft photon 
spectral index ro.i_2.4kcV and the Balmer line width FWHM(H/3) 
(Puchnarewicz et al. 1992; Boiler et al. 1996; Brandt et al. 1997; 
Williams et al. 2004). Recently Zhou et al. (2006) have reexam- 
ined this relationship using larger sample of objects and concluded 
that the anticorrelation in NLSls extends only to FWHM ~ 1000 
kms -1 , and at smaller line widths the trend appears to reverse. 

Interpretation of observational data suggests that NLSls are 
Seyfert Is, which have relatively low mass for a given luminosity 
and therefore radiate near the Eddington limit (Pounds et al. 1995; 
Hayashida 2000; Shrader & Titarchuk 2003; Botte et al. 2004); 
even super-Eddington accretion rates were suggested (Wang et al. 
1996; Collin & Kawaguchi 2004; Zhang & Wang 2006). Estimated 
Eddington ratios as well as the spectral properties may indicate 
that NLSls represent the supermassive BH analogue of Galactic 
BH (GBH) in their high states (Pounds et al. 1995; McHardy et al. 
2004; Gierlihski & Done 2006; Sobolewska & Done 2007). How- 
ever, more research should be done to confirm — or reject — this 
view, and thus the black hole mass determination is one of the key 
issues. 

In the present paper we determine the X-ray excess variance 
for a larger sample of 21 NLS1 galaxies and discuss the scaling 
properties in the combined sample of NLS1 and BLS1 objects. 

The methods of mass determination are given in Sect. [2] the 
obtained results are shown in Sect. [3] and we discuss the results in 
Sect.H 



2 METHOD 

We select a sample of NLS 1 galaxies based on the requirement that 
an object was well studied both in optical and X-ray band. Majority 



Table 2. The BLS1 sources in our sample taken from Nikolajuk et al. 2006. 



Name 


FWHM(H/3) 


ret". 


ro.l-=-2.4keV 


ref. 


3C 120 


2360 ±170 


1 






3C 390.3 


9630 ± 804 


2 


1.87 ± 0.45 


3 


Ark 120 


5536 ± 297 


2 


2.63 ± 0.29 


3 


IC 4329A 


5620 ± 200 


1 


1.71 ± 0.87 


3 


Mrk 509 


3430 ± 240 


1 


2.34 ± 0.16 


3 


NGC 3227 


5138 ± 787 


2 






NGC 3516 


3353 ± 310 


1 






NGC 3783 


3570 ± 190 


1 


3.16 ±0.63 


3 


NGC 4151 


4248 ± 516 


2 






NGC 4593 


5320 ± 610 


1 


2.49 ± 0.42 


3 


NGC 5548 


5830 ± 230 


1 


2.21 ±0.15 


3 


NGC 7469 


2650 ± 220 


1 


2.52 ±0.24 


3 


F9 


6270 ± 290 


1 


2.21 ±0.19 


3 



Column (1) list the object name. Columns (2) and (4) list the values of 
the FWHM of H/3 (broad component) in kms -1 and the values of the soft 
photon index, ro.1_2.4keV> respectively. The numbers in Columns (3) and 
(5) correspond to the following references: 1 - Nandra (2006), 2 - Peterson 
et al. (2004), 3 - Walter & Fink (1993). 

of the classical NLS1 sources have the low redshift values (z < 
0.1) and low X-ray luminosities in 2-10 keV band. Additionally 
we add two BLSls - dwarf active galaxies which are close to being 
Sy2. First object is the least luminous AGN known - NGC 4395, 
and the second object is also a low luminosity source, with jet, and 
with maser emission - NGC 4258. TableQ]shows our sample of 23 
selected sources. Types, redshifts and coordinates of the galaxies 
were taken from a catalogue of quasars and active nuclei (Veron- 
Cetty & Veron, 2006). 

We supplement it with the sample of BLS1 galaxies taken 
from Nikolajuk et al. 2006 (see Table[2j. 



2.1 Black hole mass from method based on optical 
observations 

The estimation of the black hole masses is difficult. Several meth- 
ods were developed for that purpose. We tried to collect the black 
hole mass determinations for our sample in possibly consistent way. 

The most precise method to determine A/bh is based on water 
maser emission and the black hole mass from this method was used 
for NGC 4258. 

Most of the sources were subject of AGN monitoring cam- 
paign so the reverberation method could have been used (e.g. 
Blandford & McKee 1982; Peterson & Wandel 1999, 2000; Vester- 
gaard 2002). This technique assumes that the motion of the gas 
moving around the black hole is dominated by the gravitational 
force and therefore the gaseous clouds in Broad Line Region (BLR) 
are virialized. Hence the central black hole mass can be estimated 
using the BLR radius, J?blr, and velocity of the BLR gas, «blr: 
Mbh = ^blr^blr/ G, where G is the gravitational constant. The 
radius 7?blr can be estimated from observations of a light travel 
time r. One can assume that an optical/UV continuum, which is 
produced in an accretion disc, photoionises BLR clouds and leads 
to formation of emission lines (e.g. H/3). The continuum varies in 
time and therefore the lines in BLR also vary but they respond with 
a delay r. The value of r can be thus interpreted as a light travel 
distance Rblh- Such an approach to estimations of Rblh distance 
is accurate, and it is later referred as 'reverberation method'. How- 
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Table 1. The NLS1 sources in our sample supplemented with two BLS1 objects. 



Name RA(J2000.0) DEC(J2000.0) z Type FWHM(H/3) ref. r .i^ 2 .4keV ref - 



Mrk 335 


00 h 


06 m 


19.5 3 


+20' 


12 


10" 


0.026 


NLS1 


1710 ± 140 


2 


3.10 


± 


0.05 


1 


IZw 1 


00 


53 


34.9 


+12 


41 


36 


0.061 


NLS1/QSO 


1240 


1 


3.09 


± 


0.16 


1 


TON SI 80 


00 


57 


19.9 


-22 


22 


59 


0.062 


NLS1/QSO 


970 ± 100 


2 


3.04 


± 


0.01 


1 


PHL 1092 


01 


39 


55.8 


+06 


19 


21 


0.396 


NLS1/QSO 


1790 


1 


4.3 


± 


0.3 


1 


1H 0707-495 


07 


08 


41.5 


-49 


33 


06 


0.041 


NLS1 


1050 


1 


2.25 


± 


0.25 


1 


Mrk 110 


09 


25 


12.9 


+52 


17 


11 


0.036 


NLS1 


1760 ± 50 


2 


2.17 


± 


0.27 


9 


Mrk 142 


10 


25 


31.3 


+51 


10 


35 


0.045 


NLS1 


1620 ± 120 


2 


3.15 


± 


0.12 


1 


Mrk 42 


11 


53 


41.8 


+46 


12 


13 


0.025 


NLS1 


670-865 


3.4 


2.6 


± 


0.2 


8 


NGC4051 


12 


03 


09.6 


+44 


31 


53 


0.002 


NLS1 


1170+ 100 


2 


2.84 


± 


0.04 


1 


NGC 4395 


12 


13 


48.9 


+33 


42 


48 


0.001 


BLS1/S1.8 


1500 


5 




A 






PG 1211+143 


12 


14 


17.7 


+14 


03 


13 


0.081 


NLS1/QSO 


1900 ± 150 


2 


3.03 


± 


0.15 


1 


Mrk 766 


12 


18 


26.5 


+29 


48 


16 


0.012 


NLS1 


1100 ± 200 


2 


2.75 


± 


0.13 


9 


NGC 4258 


12 


18 


57.5 


+47 


18 


14 


0.001 


BLS1/S1.9 


1200 ± 400 * 


6 




□ 






PG 1244+026 


12 


46 


35.2 


+02 


22 


09 


0.048 


NLS1 


830 ± 50 


2 


3.26 


± 


0.14 


1 


MCG -6-30-15 


13 


35 


53.8 


-34 


17 


11 


0.008 


NLS1 


1700 ± 170 


7 


2.33 


± 


0.23 


9 


PG 1404+226 


14 


06 


21.8 


+22 


23 


16 


0.098 


NLS1 


880-950 


1.4 


4.3 


± 


0.3 


1 


NGC 5506 


14 


13 


14.8 


-03 


12 


27 


0.006 


NLS1/S2.0 














Mrk 478 


14 


42 


07.4 


+35 


26 


23 


0.079 


NLS1/QSO 


1630 ± 150 


2 


3.06 


± 


0.03 


1 


HB89 1557+272 


15 


59 


22.2 


+27 


03 


39 


0.065 


NLS1 


1410 


8 


1.3 


± 


0.6 


3 


IRAS 17020+4544 


17 


03 


30.4 


+45 


10 


17 


0.060 


NLS1 


1040 


1 


2.37 


± 


0.19 


1 


Mrk 507 


17 


48 


38.4 


+68 


12 


16 


0.056 


NLS1 


960-1565 


8.4 


1.68 


± 


0.16 


1 


IC 5063 


20 


52 


02.3 


-57 


04 


08 


0.011 


NLS1/S2.0 














Akn564 


22 


42 


39.3 


+29 


43 


31 


0.025 


NLS1 


720-950 


8.1 


3.47 


± 


0.07 


1 



Column (l)-(4) list the object name, coordinates and redshifts. Column (5) list types of the galaxies. NLS1 or BLS1 mean the overall type of the galaxy. QSO 
(in the second part, after /) denotes that this source is brighter in B colour than absolute magnitude -23 mag (Mb < —23). In other case of remaining NLSls 
or BLSls we have galaxies that are fainter than -23 mag. Columns (6) and (8) list the values of the FWHM of H/3 (broad component) in kms -1 and the 
values of the soft photon index, ro.i-2.4keV. respectively. Both values were taken from literature. For NGC 4258, the sign * denotes FWHM of polarized 

H/3. The sign A in Column (8) denotes that T (above 1 keV) varies significantly on time-scales of 1 year or less (r ~ 0.6-1.72; Moran et al. 2005). 
The □ indicates that the nuclear component is visible above 2 keV only (Fiore et al. 2001). The numbers in Column (7) and (9) correspond to the following 
references: 1 - Leighly (1999b), 2 - Grupe et al. (2004), 3 - Bian & Zhao (2003), 4 - Wang & Lu (2001), 5 - Kraemer et al. (1999), 6 - Barth et al. (1999), 7 - 

Nandra (2006), 8 - Boiler et al. (1996), 9 - Walter & Fink (1993). 



ever, it requires a long observational campaign, which are time con- 
suming. In order to save that time, the radius .Rblr can be es- 
timated from the empirical relationship between the size and the 
monochromatic luminosity at 5100A (Kaspi et al. 2000): 



Rblr 



(32.9t?;°; 



AL A (5100A) 



10 44 ergs" 



o.7()()±o.(m 



light — days . (1) 



Such approach is called 'scaling method' and it was used by many 
scientist (Wandel 1999; Kaspi et al. 2000, 2005; Bentz et al. 2006; 
McGill et al. 2008). As to the velocity dispersion, it can be ex- 
pressed as «blr = /xFWHM, where FWHM is the width of 
e.g. H/3 line and the geometrical factor / depends on the geom- 
etry and the details of the gas kinematics in BLR (Krolik 2001; 
McLure & Dunlop 2001; Collin et al. 2006). Many authors assume 
the isotropy of the gas motions in BLR. In such case the factor 
/ is equal to y/3/2 (Netzer 1990, Kaspi et al. 2000) and the ve- 
locity can be written as d B lr = (V3/2)FWHM(H/3). Onken et 
al. (2004) determined new value of the factor /. They have used 
the correlation between the black hole mass and the stellar veloc- 
ity dispersion in bulge mainly in BLS1 galaxies. They found out 
/ ~ V5.5/2 = 1.17. Such value of / was adopted by Peterson et 
al. (2004) and Peterson et al. (2005) in the estimation of BH masses 
in the reverberation method. We adopt those mass measurements 
when available (i.e. for sources Mrk 335, Mrk 110, NGC 4051, 
NGC 4395 and PG1212+143). Several other objects (I Zwl, TON 
S180, PHL 1092, 1H 0707, Mrk 142, Mrk 42, Mrk 766, PG 1244, 
MCG -6-30-15, PG 1404, Mrk 478, IRAS 17020, Mrk 507, Ark 



564) had mass determination based on the old / coefficient used 
by Kaspi and those masses were corrected for the two coefficients 
ratio, i.e. by the factor |^ = 1-83 in order to obtain Mbh.l val- 
ues. The original values and the rescaled values are both given in 
Table [3] When we have a choice between the black hole masses 
from Wang & Lu (2001) or Bian & Zhao (2003) (both masses es- 
timated through the scaling method based on relationship between 
BLR radius and optical luminosity at 5100 A) we have decided to 
take A/bh, l from the later paper. We believe that masses estimated 
in 2003 are determined more precisely that in 2001. 

For two objects (NGC 5506 and in IC 5063) the black hole 
mass was estimated from stellar velocity dispersion, and the deter- 
mination of the mass for HB89 1557+272 is outlined below. 

All mass values and the references are collected in Table[3] 



2.1.1 The case of HB89 1557+272 

No black hole mass estimate for HB89 1557+272 can be found in 
the literature, and the object was not monitored. Therefore, we ob- 
tain a simple order of magnitude estimate in a following way. From 
Malkan, Margon & Chanan (1984) we take the total bulge mag- 
nitude in r (6650A; defined by e.g. Thuan & Gunn 1976) colour. 
This value m r , bulge is equal to 17.05. Then we convert this lu- 
minosity to the bulge magnitude in B (4400A) colour. For this 
purpose, we assume that the bulge spectrum in HB89 1557+272 
is the same as in M31. Using it (see Wamsteker et al. 1990) we 
obtain the ratio F 1/ (6650A)/F 1/ (4400A) = 1.24. We convert 
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Table 3. The black hole masses taken from literature and/or corrected by us 
for NLS1 galaxies and two BLS1 objects. 



Name 



M old 
(TO 6 M ) 



■Mbh, l 
(10 6 M ) 



ref./ 
meth. 



Mrk 335 




14.2 ±3.7 


1/r 


IZw 1 


18.20t!;j$ 


33-4±?; 6 


2/s 


TON SI 80 


11.6-12.3 


21.3-22.5 


3,5/s 


PHL 1092 


160 ± 60 


293 ± 110 


4/s 


1H 0707+495 


2.04 


3.74 


5/s 


Mrk 110 




25.1 ±6.1 


1/r 


Mrk 142 


4.7 


8.62 


3/s 


Mrk 42 


0.4-10.0 


0.73-18.3 


5,3/s 


NGC4051 




1.91 ± 0.78 


1/r 


NGC 4395 




0.36±g;li 


6/r 


PG 1211+143 




146 ± 44 


1/r 


Mrk 766 


0.83-4.3 


1.52-7.88 


5,3/s 






Qfl _L 1 


//111 


PG 1244+026 


1.3 


2.38 


3/s 


MCG -6-30-15 


1.55 ±0.30 


2.84 ±0.55 


5/s 


PG 1404+226 


6.6-10.0 


12.1-18.3 


5,3/s 


NGC 5506 




88 


8/d 


Mrk 478 


18.7-21.9 


34.3-40.2 


3,5/s 


HB89 1557+272 




4.17 




IRAS 17020+45 


5.9 


10.8 


3/s 


Mrk 507 


11.6 


21.3 


3/s 


IC 5063 




55 


9/d 


Ark 564 


1.2-2.9 


2.20-5.32 


5,3/s 



Column (1) lists the object name. Column (2) shows the masses of black 
holes obtained by using the scaling method (based on relationship between 
BLR radius and optical luminosity). Those values M° ld have been cor- 
rected by us and shown in Column (3) (see text). Column (3) also lists 
A/bh values of those sources for which masses were obtained from other 
methods (see the letter in Column (4)). In case of HB89 1557+272 their 
M.;.{ was calculated (see text for details). Values Mbh, l were taken for 
further analysis. The numbers in Column (4) correspond to the following 
references: 1 - Peterson et al. (2004), 2 - Vestergaard (2002), 3 - Wang & 
Lu (2001), 4 - Dasgupta et al. (2004), 5 - Bian & Zhao (2003), 6 - Peterson 
et al. (2005), 7 - Herrnstein et al. (1998), 8 - Papadakis (2004) 9 - Woo & 
Urry (2002). The letter in Column (4) shows the method of estimations of 
Mbh: r' - the reverberation mapping technique, V - the scaling method, 
'd' - stellar velocity dispersion, 'm' - maser). 



""ir.buige mt0 m B, bulge based on this information and formulas 
from Zombeck (1990). We obtain that ma, bulge = 17.73 mag. In 
the third step, we calculate Lb, bulge- Due to small redshift of this 
object (z=0.064625) we assume that in our universe the cosmolog- 
ical constant A = 0. The Hubble constant and the deceleration 
parameter are Ho = 75 km s -1 Mpc -1 and q = 0.5, respectively. 
Under those conditions we obtain the bulge luminosity in the B 
colour Lb, bulge = 1.28 x 10 43 erg s" 1 . Kormedy & Gebhardt 
(2001) have presented the correlation between A/bh and LB.buige: 



M B 



0.78 x 10" 



3, bulge 



10 10 L B 



(2) 



where Lb,q denotes the Sun luminosity in B bancfl Hence the 
black hole mass obtained from the formula is equal to 5.00 x 
10 7 AfQ. Kormedy & Gebhardt's relationship is worth only for 
normal galaxies and BLS1. Veron-Cetty & Veron (2006) classify 



HB89 1557+272 as ordinary NLS1 sources, because it luminos- 
ity in B band is low. However, Malkan et al. (1984), Hutching et 
al. (1984) or NEE0 classify this object as quasar with the narrow 
emission lines. We assume the narrow line type of HB89 1557+272 
and have to correct obtained Mbh in order to estimate true value 
of it. Mathur, Kuraszkiewicz & Czerny (2001) have mentioned that 
the BH to bulge mass ratio for NLS 1 galaxies is of order of 0.00005 
and 0.0005 for NL quasars. For Wandel's (1999) sample of galax- 
ies this ratio M B H/M b uigc = 0.0003 for BLS1 and 0.006 for normal 
galaxies and quasars. We adopt those values in our analysis. If we 
assume that (i) the bulge masses in NLS1 and BLS1 galaxies are 
the same (M^f 1 = M^^), (ii) the BH masses in NLS Is are si- 
multaneously smaller than in BLSls (A/bh S1 < A/|h S1 ), and (iii) 
there is only a systematic shift in equation l[2} between those two 
kinds of Seyfert 1 galaxies, we obtain that ALJh 31 = 6 x A/bh S1 



and M 



BLSl/QSO 
BH 



= 12 x M 



NLS1/QSO 
BH 



Therefore, we should di- 



vide the value 5.00 x 10' by 6 in the case if HB 89 1557+272 would 
be NLS 1 and by 12 if the galaxy would be NL quasars. We choose 
the last one factor and the BH mass in HB89 1557+272 is A/bh= 
4.17 x W 6 Mq. We adopted this value in our further analysis (Ta- 
ble [3). We must note here that the value of Mbh calculated by us 
has a large error and may be underestimated. 



2.2 Determination of the X-ray excess variance 

We determine the X-ray excess variance for the sources in our sam- 
ple using the approach which was developed by Nikolajuk et al. 
(2004) with the aim to use for mass determination in BLS 1 galax- 
ies. The approach is based on the use of only the high frequency 
tail of the power spectrum and on the specific way of combining 
results for several lightcurves for a given source. 

Each AGN emit a variable X-ray radiation observed in the 
hard (2-10 keV) band. From such lightcurve we can calculate the 
normalized excess variance which is defined as (e.g. Nandra et al. 
1997; Turner et al. 1999): 



1 N 

— [fa - S) 2 - (Sxif] 



(3) 



Here, N is the number of data points, Xi and 5xi are the flux and its 
error, respectively, x is the unweighted, arithmetic mean of Xi. cr„ xs 
is in units of (rms/mean) 2 . The important value is not the excess 
variance itself but actually the ratio of the excess variance to the 
duration of the lightcurve since the power spectrum has the slope 
~ 2 at high frequencies. This is clearly seen from the formula used 
to calculate A4bh for BLS1 galaxies: 

-C T -^, (4) 



A/ : 



where T is the duration of the X-ray lightcurve (in seconds), At 
is its bin size (in seconds) and C is the constant. The requirement 
of the method is that the excess variance is measured at the high 
frequency tail of the power spectrum, i.e. that the length of a single 
lightcurve, T, is significantly shorter than the inverse of the high 
frequency break of the power spectrum, i-v (i.e. T < —). The 
same condition is important for any scaling relations so we follow 
it here. We would like to emphasise here that we need not to know 
exact value of T. We chose T based on information that our T 
should be shorter than break time, Tb r . This time, which is simply 
equal to inversion of i^, r (i.e. Tt, r = l/^ r ), is shown in Table [4] 



1.93 x 10 33 i 



NASA/IPAC Extragalactic Database 
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In case of those NLS 1 sources, for which we do not have any in- 
formation about Tb r , we estimate this time in following way. We 
use the relationship between the black hole mass and the break fre- 
quency in BLS1 galaxies (it, r = 15/(A/bh/Mq) Hz) derived by 
Papadakis (2004). In our calculations we assume Mbh = Mbh,l- 
The break frequencies in NLSls are higher by a factor 10-30 than 
in BLS1 sources (e.g. Papadakis 2004). Therefore, we multiply 
Papadakis' relationship by a factor D — 10 — 30, i.e. our esti- 
mated Tb r = A'/bh, l/15D s. Generally, the value D is equal to 
30. Sometimes the calculated excess variance had been less than 
zero and in those cases we decreased D in order to obtain longer 
T. We have decreased D from 30, unless the values of ct^xs become 
positive, but still D > 10. We choose now T = [0.7 - 0.95] xT br . 

More accurate mass determination can be obtained if more 
measurements are available for a given source. If a single but long 
observation is available, it can be split into several parts. Instead 
of direct computation of an average excess variance we use a pro- 
cedure which allows for different duration of the lightcurve parts 
and higher determination accuracy. The separate measurements are 
combined in the following way: (i) we obtain the kth individual 
ratios: 



Table 4. X-ray observations and selected lightcurves for NLS 1 sample and 
two BLS1 objects. 



A k = 



(°"nxs)fc 

(T-2At) k 



(5) 



M- 



(7) 



(ii) we fit this set of the individual Ak and calculate the mean value 
of (A) as weighted mean from the minimum chi-square fit, in order 
to estimate the final value of the X-ray excess variance 

*L s = {A){T-2At) , (6) 

where (T — 2 At) is the mean duration of a single lightcurve. The 
value of (A) was also used to calculate the black hole mass from 
the excess variance (Nikolajuk et al. 2006) 

C 

' ~ (A) ' 

Gierlinski et al. (2008) have shown that the method applies very 
well to galactic sources at their hard states and to Seyfert 1 galaxies, 
and the proportionality constant, C, does not depend on the source 
luminosity state. 

The value of the scaling constant was determined to be C = 
0.96 ± 0.02 MqS" 1 (Nikolajuk et al. 2004), if based on the black 
hole mass of Cyg X-l equal to 10 Mq (e.g. Orosz 2003, Herrero 
et al. 1995; Nowak et al. 1999; Gierlinski et al. 1999). Nikolajuk et 
al. (2006) obtained 1.92 ± 0.5 Mq s" 1 , again based on Cyg X-l 
scaling but taking the Cyg X- 1 mass of 20 ± 5 Mq after Ziolkowski 
(2005). Recent study of several galactic sources by Gierlinski et al. 
(2008) gave the value 1.24 ± O.OfJf]. In the present paper we adopt 
the value of 1.92 from Nikolajuk et al. (2006) for better consistency 
with the previous study of BLS 1 galaxies. 

As it was mentioned by Nikolajuk et al. (2006), the accuracy 
of the excess variance and M BH 2 determination depend on the 
effect of statistical error of variance measurement and power leak- 
ing from long time-scales. The errors are estimated by performing 
Monte Carlo simulations for each source separately. The simulation 
generated a few hundred sets of artificial data. From each lightcurve 
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T br (ref.) 


Proposal 


Mean 
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name 


[ks] 


details 


[cnts/s] 
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4.1 19 


20 






R/P40155 
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NGC 5506 
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Mrk 478 
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5 


HB89 1557+272 




A/81004000 
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0.559 
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IRAS 17020+45 




A/73047000 
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11 


Mrk 507 
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0.042 
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IC 5063 
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A/71030010 


0.124 
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Ark 564 
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R/P10291 


2.035 


19 



Column (1) lists the object name. Column (2) shows the power spectral 
density (PSD) break time-scales, which are connected with the frequency 
break (Tb r =l/^b r )- The number in parentheses correspond to the following 
references: (1) McHardy et al. (2004); (2) Vaughan et al. (2005); (3) 
Vaughan & Fabian (2003); (4) Markowitz & Uttley (2005); (5) McHardy et 
al. (2005); (6) Uttley & McHardy (2005); (7) Papadakis et al. (2002). Col- 
umn (3) shows the observation details. The first letter refers to the satellite 
name (A - ASCA, R - RXTE), the number refers to proposal number. In 
Column (4) are shown the mean fluxes observed in the lightcurve (in count 
per second units). Column (5) lists the number of short lightcurves, Ni c , 
subtracted from whole proposal lightcurve. N\ c also means the number of 
the values (cr^ xs )fc, which were used in order to estimate A7 BH CT 2. 



the least squares method, the (^4). Finally, the 90% error of {A) 
was calculated due to having a distribution of the hundreds \values 
of them. The mass error was determined from the error of (A). 



3 Note the difference in the definition of the constant C between the present 
paper and Gierlinski et al. which leads to a difference in units but not in the 
numerical value: Cni^o = Coier x v o> where vq = 1 Hz. Hence, the 
numerical value of the coefficients are the same in units Mq s — 1 and M© 
Hz" 1 . 



3 RESULTS 

Since NLS 1 galaxies are known to be relatively more variable than 
BLS1 galaxies (e.g. Hayashida 2000, Czerny et al. 2001; McHardy 
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Table 5. The mean time of observations and the mean normalized excess 
variances the whole sample of galaxies. 



Table 6. The black hole masses. 



Name 



(T) 
[si 



T 2 

n xs 



(rms/mean) 



Mrk 335 
IZw 1 
TON SI 80 
PHL 1092 
1H 0707+495 
Mrk 110 
Mrk 142 
Mrk 42 
NGC4051 
NGC 4395 
PG 1211+143 
Mrk 766 
NGC 4258 
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IC 5063 
Ark 564 
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x 10~ 3 
x 10~ 3 
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x 10~ 3 
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x 10~ 3 
x 10" 3 
x 10~ 3 
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IC 4329A 
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NGC 3227 
NGC 3516 
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NGC 4151 
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NGC 5548 
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x 10 
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x 10 
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F303 

PG 0026+129 
Mrk 586 
Mrk 1040 
NGC 985 
Mrk 279 
Mrk 841 
Mrk 290 



39936 
39936 
39936 
39936 
39936 
39936 
39936 
39936 



(6.72 ±6.03) 
(1.31 ± 1.92) 
(2.57 ±0.75) 
(1.20 ±0.65) 
(3.47 ±1.76) 
(2.32 ±0.84) 
(1.14 ±0.93) 
(4.11 ±2.15) 



x 10~ 3 
x 10" 3 
x 10" 2 
x 10" 2 
x 10~ 3 
x 10 -3 
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43.07 (a 
43.35 (a 
43.58 (a 
44.15 (b 

42.49 (b 
43.80 (c 
43.17 (a 
42.12 (b 
41.21 (a 
39.99 (a 
43.63 (b 
42.73 (a 
40.52 (d 
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42.72 (a 
43.03 (b 
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43.50 (a 



43.62 (b 
42.87 (e 
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43.95 (d 
44.00 (d 
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43.59 (a 
44.03 (a 
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43.03 (a 

44.53 (a 
44.07 (a 
42.40 (a 
43.50 (a 
43.66 (a 

43.54 (a 
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First part of this table concerns the values of sources considered in this pa- 
per. Second part shows the values obtained for BLS 1 galaxies (sample taken 
from Nikolajuk et al. 2006). Galaxies and their appropriate values shown in 
third part are taken from O'Neill et al. 2005. Those all objects and sev- 
eral galaxies from part one and two have been considered by Liu & Zhang 
(2008) in their analysis. Column (1) lists the object name. Column (2) shows 
the arithmetic mean time of observations. Col (3) gives the normalized ex- 
cess variance with errors. Those values have been calculated taking mean 
(T — 2 At) and (A) as the product of those two quantities (see equation|6j. 
The errors have been obtained from Monte Carlo simulations of (A) (see 
text). Column (4) shows logarithm of the luminosity in the 2-10 keV band. 
The cosmology is Ho = 75 kms _1 Mpc _1 and <jo = 0.5. References in 
this column: (a) - O'Neill et al. 2005, (b) - Vaughan et al. 1999, (c) - George 
et al. (2000), (d) - Merloni et al. 2003, (e) - Bassani et al. (1999). 
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Mrk 507 
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Ark 564 
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Column (1) lists the object name. Column (2) shows the masses of black 
holes without errors taken from Table|3] Column (3) shows the black holes 
masses obtained from the variance method. The ratios of the two masses are 
given in Column (4). 



et al. 2004; Papadakis 2004), we do not expect the standard X-ray 
excess variance method of mass measurement, with the same value 
of the constant C, to work for these objects. 

Instead, we search here for the discrepancies between variabil- 
ity based method and other methods of mass measurements, and we 
analyse the dependence of these discrepancies on the source prop- 
erties. 

The details on the lightcurves analysed for the purpose of this 
work are given in Tableland the values of the X-ray excess vari- 
ance for a whole sample of both NLS 1 and BLS 1 galaxies are given 
in Tabled 

We first check single parameter relations between the X-ray 
excess variance and the black hole mass (taken from optical de- 
termination) and the X-ray luminosity. We renormalize the excess 
variance to the value expected for the duration of the lightcurve 
equal to 40 000 s, i.e. cr| = cr^s X (40000/ (T)). The results are 
plotted in Figs.[T]and|2] 

There is much better correlation with mass than with lumi- 
nosity in the whole combined sample of NLS1 and BLS1 galaxies 
(correlation coefficient equal to -0.69 and -0.39 correspondingly). 
However, the correlation with the luminosity is basically driven by 
the BLS 1 objects. If only NLS 1 are selected, the correlation of the 
excess variance with the black hole mass is stronger than with the 
X-ray luminosity (-0.44 and -0.27 correspondingly). 

The correlation with the mass suggest that the X-ray excess 
variance can possibly be used for mass determination in NLS1 
galaxies. We therefore calculate the masses of the black holes in 
our sample using the equation Q. 

© 0000 RAS, MNRAS 000, 000-000 




Figure 1. The renormalized excess variance against the X-ray luminos- 
ity in the 2-10 keV range. The renormalized excess variance cr| = 
fnxs x (40000/ (T)) i.e. to the value expected for the duration of the 
lightcurve equal to 40 000 s. The whole sample of galaxies is presented 
in TablefS] We enlarge our sample by a few galaxies taken from O'Neill et 
al. (2005) in order to be consistent with sample of Liu & Zhang (2008). 



Figure 3. Black hole masses estimated from X-ray variance 
method,Af BH CT 2, versus masses taken from literature, A/bh.l- The con- 
tinuous line shows the relation A/bh l = M BH CT 2 ■ The dot line shows the 
relation Mbh, l = 20 X M BH a 2. We can see that small subset of NLS1 
objects which do not require any rescaling like BLS1 galaxies. The arrow 
close to NGC 4395 indicates that Mbh. l of those galaxy may be smaller 
in agreement with Filippenko & Ho (2003) (A/ BH , L ~ 10 4 - 1O 5 M ). 
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Figure 2. The same as in Fig.[TJbut versus the black hole masses taken from 
optical determination. 



The values of the black hole masses derived from the standard 
excess variance, M BHcT 2, the values of the masses based on other 
methods (reverberation, stellar velocity dispersion, etc), Mbh,l> 
and the ratio of the two values for each source in our sample of 
NLS1 galaxies, r = ,. BH|L , are given in Table|6] 

M BH,cr2 

In Fig. [3] we show the plot of the mass from variability vs. 
mass from the literature for our sample of galaxies and for BLS1 
galaxies from Nikolajuk et al. (2006). The BLS1 sources cluster 
along the continuous line marking the ratio 1. Apart from two 
sources, the value of r for all objects is within the limits 0.3 - 3. 
The source with the smallest mass, NGC 4395 is shifted from the 
line (r = 5.81 for this source). We have to note that our obtained 
value (6.2 x 1O 4 M0) is similar to other values obtained from other 
methods based on X-ray variability: ~ 6 x 10 4 Mq (Vaughan et al. 
2005; M BH - P{y) xv= 10~ 3 method); 1O 4 -1O 5 M (Shih et al. 
2003; A/bh — method). Other black hole mass determinations 
for this source also gave lower values than the direct reverberation 
(photoionization, Kraemer et al. 1999; -Rblr scaling with luminos- 
ity, Filippenko & Ho 2003) have mentioned Af BH in NGC 4395 of 
order of ~ 10 4 -10 5 Mq. In such case the ratio r would be equal 
to 0.6-1.61, respectively. On the other hand, the upper limit from 
stellar dispersion velocity in the central 3.9 pc (6.2 x 10 6 A/q, Fil- 
ippenko & Ho 2003) is still consistent with the reverberation value. 
In the case of NGC 4258 the BH mass is estimated precisely from 
maser emission method. Our determination gives the upper limit of 
A/bh 15.2 x 10 6 M Q and hence the lower limit of r = 2.5, so the 
agreement is not quite satisfactory. 

The NLS1 sources show larger dispersion than BLS1 sources. 
Most of them cluster along the dotted line, marking the ratio of 20. 
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Figure 4. Two histograms showing the number of our sources (NLS1 + 
BLS1) in a particular bin of logarithm of the mass ratio Mbh l/^bh <t 2 - 
Histogram, which is drawn by long dash line, concerns all sources in our 
sample. Shaded histogram, drawn by continuous line concerns objects with 
masses Mbh. l estimated based on the reverberation method or by maser 
method (only NGC 4258). We can see bimodality in both cases. 



Figure 5. The Full Width at Half Maximum of H/3 line produced in Broad 
Line Region in NLS1 galaxies versus black hole masses Mbh, l- The con- 
tinuous, dot-short-dash and dot lines represent accretion rate in Eddington 
units, m (respectively 1, 0.1 and 0.01). 



This simply reflects the fact that in both plots of o"n XS (Figs. [T] and 
[2]l the NLS 1 objects lie on average higher, as expected. 

The existence of two characteristic values of the ratio is al- 
ready suggested by Fig. [3] but it is seen more clearly from the di- 
rect plot of the mass ratio. In Fig. [4] we present the results on the 
variability level in our sample in a form of two histograms show- 
ing the number of sources in a particular bin of the A/bh. l and 
M B h,o-2 mass ratio. First histogram (drawn by long-dashed line) 
considers all objects in our sample of NLS1 and BLS1 galaxies. 
Second histogram (drawn by continuous line and shaded) shows 
galaxies which A/bh,l are taken from the reverberation or the 
maser method (for NGC 4258) only. Both histograms are simi- 
lar. The distributions show two separate peaks, one close to unity, 
and one close to the factor 20. BLS1 galaxies occupy the first peak 
while most of NLS 1 galaxies make the second one. 

Our data come from two different instruments (ASCA and 
Rossi-XTE; see Table |4]l which cover slightly different energy 
range we checked whether this may have a strong impact onto our 
results. For MCG -6-30-15 we had both ASCA and Rossi-XTE 
data, so we calculated the X-ray excess variance for this source 
independently from both instruments. The ASCA variance was by 
a factor 1.957 higher than Rossi-XTE variance which is consis- 
tent with ASCA being sensitive in softer energy band and with 
the general trend of lower normalization of the power spectrum in 
higher energies (e.g. Uttley & McHardy 2005). We then made an 
experiment and multiplied all variance values from Rossi-XTE by 
this factor and replotted the histogram. It still did not change and 
showed apparent bimodality. 

However, there exists a small subset of our NLS1 sources 
which does not require that multiplication and the value of r is the 



same like in BLS1 cases (MCG -6-30-15, Mrk 42, Mrk 507, Mrk 
766 and quasar HB89 1557+272). 

Since the behaviour of this small subset of NLS1 is puzzling 
and may help to understand the origin of the variability we analyse 
the possible cause of this behaviour. 

The location of these slowly variable NLS1 on the 
FWHM(H/3) versus the mass diagram is shown in Fig. [5] The er- 
rors are large but the sources do not seem to be located differ- 
ently from the other NLS1 objects. In Fig. [5] we also indicate the 
lines corresponding to the fixed dimensionless accretion rate, or 
the Eddington ratio, ?h = M/AfEdd. We have calculated m using 
equation taken from Czerny, Rozanska & Kuraszkiewicz (2004): 
FWHM = AM^'rh' - 35 , where A is constant. They derived 
this formula using Kaspi et al. (2000) relationship for .Rblr and 
assuming the geometrical factor / = y3/2. We have corrected 
Czerny et al. formula by taking new estimated by Bentz et al. 
(2007) relationship Rblr oc [AL A (5100A)] ' 54 and assuming 
/ = V5.5/2. Therefore, the corrected formula using by us in this 
paper is: 



FWHM(H/3) 



1650 kms- 



7 Af BH 
V10 8 M G 



(8) 



In Fig. [5] the division into BLS1 and NLS1 objects corresponds 
roughly to the division into m < 0.1 and m > 0.1 for a lO 7 A/0 
black hole. For other values black hole masses the exact classifica- 
tion according to the spectra (NLS1 versus BLS1) and according 
to the Eddington ratio do not overlap. The issue of the usefulness 
of the current definition of NLS 1 was already raised by Mathur & 
Grupe (2005a,b) in the context of comparison of optically-selected 
and X-ray selected NLS 1 . 

In Fig.[6]the NLS1 objects which have slower variability than 
most NLS Is are naturally located at the mass ratio ~ 1, but their 
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* NLSl 

. * NLSl like BLS1 

• BLS1 



l - 




HB89 1557+272 



L °g ( M BH, l/ M = 



Log rn 



Figure 6. The dimensionless accretion rate, m, plotted versus the ratio 
A/bh l/M bh ct 2. The dotted vertical lines show how many times should 
be multiplicated M BH a i in order to obtain M ;-{, l- The continuous line 
shows the best fit to our sample of galaxies. The dot-short dash line is the 
McHardy et al. fit. See Sect. |3.1| for details. 



Figure 7. The soft X-ray slope, ro.i— 2.4keV> plotted versus the dimen- 
sionless accretion rate, m. The continuous line shows fit to our sample of 
galaxies. The long-dash line shows fit to the same sample subtracted by the 
subset of 5 'NLSl like BLS1' objects. 



Eddington ratios span rather wide range. This may be related to 
relatively large errors for a single object. 

The formula used by us derives the value of the dimensionless 
accretion rate partially from the black hole mass, and relies on the 
scaling of the Broad Line region so it is not independent from the 
mass ratio in the statistical sense. Therefore, the result may be bi- 
ased by the underlying correlation. However, at present this is the 
simplest way to see the trend with the Eddington ratio. We will 
discuss Fig.[6]again in Sect. 13. 11 

In order to gain additional insight which may not be biased by 
determination of m, we consider the dependence of the variability 
on the soft X-ray slope. 

The spectral slope in 0.1-2.4 keV range in our sample gener- 
ally correlates with the dimensionless accretion rate (see Fig. |7J- 
For two slowly variable NLSl the X-ray slope is flatter than ex- 
pected from the overall correlation while for other three sources 
are close to the best fit line. The continuous line represents our 
best fit to our sample (correlation coefficient of this fit is 0.515): 
ro.i-2.4k e v = (0.48 ± 0.16) log m + (3.08 ± 0.18). The long- 
dash line shows fit to the same sample but without the subset NLS 1 
like BLS1 galaxies. This fit is r .i- 2 .4kcv = (0.53 ± 0.14)log 
m + (3.26 ± 0.16) with correlation coefficient 0.650. 

Now we present the relation between the soft X-ray slope and 
the mass ratios for our sample of NLS 1 galaxies (see Fig. [8j. A 
strong trend with a soft X-ray slope is seen again. Since the mea- 
surement of the soft X-ray slope is more direct than the estimate of 
the accretion rate or bolometric luminosity, this result indeed sup- 
ports the claim of McHardy et al. (2006) that the black hole mass 
alone does not determine the variability properties if we include 
objects with steep soft X-ray spectra into our sample. The change 
of variability properties with the change of the spectral slope is 



lix 



>0x 




0.5 1 1.5 

Log (M BH , L /M BHi ) 



Figure 8. The soft X-ray photon index, ro.i— 2.4keV. plotted versus the 
ratio A/bh. l/^bh a 2 f° r NLSl galaxies. The dotted vertical lines rep- 
resent the same as in the Fig. [6] The dashed line shows arbitrary bor- 
der ro.i_2.4kcV = 2. Horizontal arrow near HB89 indicates that the BH 
mass in this quasar may be underestimated. Vertical arrows mark that in 
agreement with Branduardi-Raymont et al. (2001) or Lee et al. (2001) the 
values of r .i_2.4k e V m MCG -6-30-15 and Mrk 766 may be lower. 
Continuous line represents the best fit ro.i_2.4keV= (0.78 ± 0.21) 
log(A/BH, l/Mbh CT 2 ) + (2.00 ± 0.27) with correlation coefficient equal 
to 0.668. 
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not surprising, since the underlying mechanism may be different as 
well. 

If indeed there are two separate mechanisms of the hot plasma 
formation, separate for soft spectra and separate for hard X-ray 
power law, then the dependence of variability level on the sources 
properties may be rather bimodal, reflecting a rapid state transition, 
instead of showing a continuous dependence on the Eddington ra- 
tio, as proposed by McHardy et al. 

Our sample may not still be large but at present the distribution 
looks rather bimodal, with sources concentrating around the ratios 
~ 1 and ~ 20 instead of continuously covering the whole diagram 
between some minimum and maximum values, as it would be ex- 
pected from a continuous dependence on rh. The corresponding 
histogram of m does not show such a bimodality. However, a his- 
togram of slopes in our sample again shows a bimodal behaviour, 
with objects clumping around r ~ 2 and r ~ 3. This supports the 
view that indeed a spectral change from hard to soft X-ray slope is 
underlying the change in variability properties of AGN. 

3.1 Multidimensional fits 

In the previous section we concentrated on the relation between the 
X-ray excess variance with the black hole mass but there were sug- 
gestions that the variability may correlate both with the mass and 
one more parameter like X-ray luminosity or bolometric luminos- 
ity. 

Since the division of objects according to the Eddington ratio 
seems to be based on better physical grounds we test whether the 
dependence of the variability level on this ratio can account for the 
issue of enhanced NLS 1 variability. Such an idea was put forward 
by McHardy et al. (2006) for a small sample of AGN. 

The dependence of the mass ratio, or the variability enhance- 
ment, on the dimensionless accretion rate for our sample is shown 
in Fig. [6] We have calculated rh using equation $8$ with the data of 
FWHM(H/3) and M BH , l taken from Tables[Tj3] Our best fit (solid 
line in figure) to the whole sample (without NGC 5506 i IC 5063 
for which we have no rh estimate) gives 

log rh = a log r + b, (9) 

with a = 0.79 ± 0.18 and b = -1.43 ± 0.14. The correlation is 
statistically significant (correlation coefficient 0.712). 

We can compare this result with the result of McHardy et al. 
(2006). If we assume that the high frequency slope of the power 
spectrum is equal 2, and the normalization of the frequency multi- 
plied power spectrum is independent from the object we can relate 
their result based on frequency break to our result involving the 
high frequency tail. With this assumption, the result of McHardy 
et al. (2006) imply we should expect a linear trend with the slope 
~ 1, or more precisely, a relation 

A — B 

log rh = — — — log r + const, (10) 
B 

where A = 2.17±g;|§ and B = 0.90±2;| are the coefficients from 
the McHardy et al. fits. Our larger sample is roughly consistent 
with this result since our slope a — 0.79 and the slope implied by 
McHardy et al. (2006) results is a =1.4. Since the dependence on 
the dimensionless accretion rate implied by the results based on the 
frequency break is even steeper it means that the overall trend of 
an increase of the mass discrepancy between the two methods with 
the Eddington ratio seems to be real although determination of this 
ratio is highly uncertain. 

We also show in this figure the fit calculated by McHardy et 



al. (2006) made for 10 objects (i.e. 5 BLS1 and 5 NLS1). This fit is 
represented by the dot-short-dash line with the slope (A — B)/ B = 
1.41 and const value chosen in arbitrary way. Close to this fit lie 
our second fit (not shown in Fig [6]>. This fit is calculated for our 
sample of the same objects as in McHardy et al., but without NGC 
5506. The value of slope a of our second fit (a = 1.07 ± 0.35) and 
correlation coefficient (0.780) indicates that our approach is indeed 
equivalent to McHardy et al. (2006) analysis. 

However, such a two-dimensional fit simply joints the extreme 
values of our mass ratios in Fig. [7] and does not provide an expla- 
nation of the bimodal distribution seen in Fig. [4] 

We also repeat the two-dimensional fit to both the black mass 
and the X-ray luminosity proposed by Liu & Zhang (2008). If we 
allow for a general relation between the X-ray excess variance in 
our sample of 43 objects (renormalized to 40 000 s lightcurves) and 
the X-ray luminosity 

l0g(Tn X a = AlogM B H,L + B log L 2 - lOkcV + C (11) 

we obtain the fit parameters A — — 0.841q'2q, B = 
— 0.0088±3;g||, C = 4.5±J;£. This sample includes the objects 
with the mass taken from the scaling relation since the number of 
sources with direct reverberation mass measurement is low. How- 
ever, it seems that the most important effect is the careful measure- 
ment of the X-ray excess variance. If we take only the excess vari- 
ance measurements obtained by O'Neill et al. (2005) and used by 
Liu & Zhang, and supplement them with our measurements only 
for sources which were not included in their study, the result is 
widely different: the dependence on the black hole mass is weak 
(A = —0.058) while the dependence on the X-ray luminosity is 
strong (B = —0.36) irrespectively of the use of scaled values of 
the black hole mass. Our measurements of the X-ray excess vari- 
ance for a number of objects differ from the values obtained by 
O'Neil et al. (2005) due to another way of averaging and different 
length of single curves. 



4 DISCUSSION 

The X-ray variability of the BLS1 galaxies at high frequencies de- 
pends only on the black hole mass, and the same is true for the 
galactic sources in their hard states (Gierlinski et al. 2008). There- 
fore, the X-ray excess variance is a viable method of black hole 
mass determination (e.g. Papadakis 2004; Nikolajuk et al. 2004; 
Awaki et al. 2005; Zhang et al. 2005; Markowitz et al. 2006; Niko- 
lajuk et al. 2006; see also O'Neill et al. 2005). 

The issue complicates when the sources with soft X-ray spec- 
tra are considered. In the case of AGN the sources with the soft 
X-ray spectra usually belong to the NLS1 galaxy class, and usu- 
ally they are more variable than the corresponding broad line ob- 
jects (Iwasawa et al. 2000; Markowitz & Edelson 2004; Vaughan et 
al. 2005). McHardy et al. (2006) suggested that this effect should 
be taken into account in a form of continuous dependence of the 
variability rate not just on the black hole mass but also on the Ed- 
dington ratio. They have found a relation which roughly represents 
the variability properties in a combined sample of NLS 1 and BLS1 
objects. 

In the present paper we studied a larger sample of 21 NLS1 
galaxies. We confirm the dependence of the variability properties 
on the accretion rate if NLS 1 galaxies are included. However, our 
study indicated that such a parametrization of variability as pro- 
posed by McHardy et al. (2006) for AGN may not fully reflect the 
complexity of the issue. 
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We have found that the essential effect is not in the dependence 
on the accretion rate but likely in the dependence on the soft X-ray 
slope. The slope is apparently related to the character of variability, 
as is seen for example in the study of the energy-dependent frac- 
tional variability by Gierlihski & Done (2006). Typical soft NLS1 
like NGC 405 1 showed much higher amplitude at lower energies, 
below 0.5 keV than at 10 keV, indicating a pivoting trend while the 
hard spectrum NLS1, MCG -6-30-15 showed the same amplitude at 
very low and very high energy, i.e. no change of the spectral slope. 
What is more, we see an indication of a bimodal behaviour con- 
nected with predominantly bimodal distribution of the soft X-ray 
slopes of AGN in our sample while the distribution of the accre- 
tion rates seems to be continuous. This explains why Gierlihski et 
al. (2008) do not find the dependence of the variability rate on lu- 
minosity in their study of the galactic sources: the sources were 
selected to represent the hard state only and the evolutionary stages 
with soft spectra were not selected for consideration. 

The fact that it is not the luminosity change but the spectra 
change, which affects the fastest variability, hints for a different, 
or additional variability mechanism when the spectrum becomes 
softer. 

We would like to emphasise the fact that our approach in 
calculating the normalized excess variance, based on computation 
weighted mean cr nxs from the minimum chi-square fit, is more re- 
liable and gives more accurate results than direct computation. In- 
deed, the way of calculating the excess variances is very important 
and it influences on obtained results. The findings depends on it 
very strongly. From this follows differences between our and Liu 
& Zhang (2008) results. When we repeated calculation of Liu & 
Zhang taking only their values of a„ xs , we obtained results similar 
to them. 

Nevertheless, we must note that we had still too small sample 
of galaxies, the data came from two different instrument (ASCA 
and Rossi-XTE), and particularly the number of sources with mass 
determination from monitoring was far from satisfactory. In order 
to resolve the fact that existence of bimodal distribution is realised 
in nature or not we need enlarge sample, especially by NLS 1 ob- 
jects which behave in the X-ray variance like BLS1 galaxies. 
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